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Abstract 

At tree level, I discuss modifying the FNAL heavy quark formalism to 
include a twisted mass term. I find that at maximal twist the so called 
KLM factor is independent of the heavy mass. 



1 Introduction and motivation 

Although only recently developed, twisted mass QCD is already proving to be 
an excellent technique for producing accurate lattice QCD results. Twisted mass 
lattice QCD [HE] calculations have been used to test and constrain chiral pertur- 
bation theory [31 HI [5], study baryons [6], heavy-light mesons [7j, flavour singlet 
mesons [S], static-light mesons [§], and the pion form factor [10J. A key theoret- 
ical advantage of the twisted mass formalism is that the action is automatically 
0(a) improved at maximal twist [21 [TT]. Twisted mass QCD has recently been 
reviewed [T2l [T3] . 

An important part of the twisted mass program is that there are power count- 
ing arguments to understand the 0(a 2 ) corrections for light quarks [Hj. These 
are based on a Symannzik analysis that is suitable for light quarks. As the mass 
of the heavy quark increases towards the mass of the charm quark and beyond, 
then 0((aM Q ) 2 ) effects may become sizable. One way to estimate heavy quark 
mass effects is to use the FNAL heavy quark mass formalism [15] . 

The full FNAL heavy quark formalism requires tuning the terms in the heavy 
quark action. For example the clover coefficient of the clover term and the co- 
efficient of the spatial Wilson term. There have been a few numerical studies of 
the required tuning p21 [T71 [181 E3 EH EI] • However, many groups have used the 
FNAL formulation to estimate heavy mass corrections to decay constants and 
masses [22l [23l El] for the standard Wilson and clover actions. The prescription 
was to use the kinetic mass and to multiply the quark fields by the KLM factor 
in equation [T] 

Z KLM = V2ne m o (1) 

where mo is defined in 

mo = log(-L-3). (2) 



and k is the standard hopping parameter used in the clover and Wilson actions. 
The corrections to equation [T]that include non-perturbative 0(a) mass corrections 
are in [25] . 

A numerical test of the KLM factors is reported by El-Khadra et al. [16j and 
a test of the FNAL formalism for the charm mass reported by Dougall et al. [26J. 
A critical comparison of the KLM factors to renormalisation factors determined 
non-perturbatively is in [27J . The aim of this paper is to find the equivalent KLM 
factor for twisted mass fermions at tree level. 

Throughout this paper I will mostly only consider tree level perturbation the- 
ory. This is in the spirit of estimating and correcting the leading 0((arriQ) n ) 
corrections to the results of twisted mass calculations. I also don't consider the 
two doublet twisted mass formalism for including non-degenerate quarks in un- 
quenched calculations [SHUSH] • The results will be useful to analyse existing n/=2 
twisted mass calculations and twisted mass calculations that use the Osterwalder- 
Seiler action [301I31] for the heavy quarks in unquenched calculations with 2 + 1 + 1 
flavours of sea quarks. 

The theoretical foundations of twisted mass QCD [U [321 E] use a mass in- 
dependent renormalisation scheme [33], but the FNAL heavy quark formalism 
uses mass dependent renormalisation factors. The issue of mass dependent ver- 
sus mass independent renormalisation schemes is reviewed by Georgi [34J and 
Kronfeld 1351 . 



2 A brief introduction to twisted mass QCD 

I first review the twisted mass quark action in the continuum. 

S F = J d^xxil^D^ + m q + ^q7 5 t 3 )x (3) 

where r 3 is the third Pauli spin matrix in flavour space, and m q and fj, q are mass 
parameters. The fields \ and x are in the twisted basis. The quarks fields can 
be transfered to what is known as the "physical basis" by the transformation: 

if) = exp(iwy 5 T 3 /2)x i/> = xex'piiuj^T' ; /2) (4) 

where tana; = fi q /m q . After the transformation in equation HJ the twisted quark 
action in the physical basis is 

S F = j d^x^D^ + ^w? q + ^ (5) 

The lattice version of the twisted mass action in equation [3J is written down 
in the standard way. However, to prepare for applying the FNAL heavy quark 
formalism, I consider a twisted version (S^f) of the action written down by Lin 
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and Christ [Tj5]. Lin and Christ [Tj5] have a table that summaries the choices of 
coefficients used by other heavy quark formulations [121 ES] • The version of the 
twisted Wilson action used in existing numerical calculations uses ( = r s =1. I 
only use arbitary value of £ and r s in section [51 The twisted QCD action in the 
twisted basis is 

Siat,F = ^ x(z)(7oA) + CliDi - -Dl - rs^Dj + m q + ir 3 ~f 5 fi q )x(x) (6) 

X 

The derivatives are defined by 

D^{x) = l -[U,{x)^{x + £) - UfrxMx - £)] (7) 

and 

Dfa(x) = U^x)^(x + p) + Ul(x)iP(x -p)- 2i(>(x) (8) 

The special choice of uo = 7r/2 is known as maximal twist, where all of the quark 
mass is in the 75T 3 term. At maximal twist there are no 0(a) corrections to 
the continuum result [2l [TTJ. Achieving maximal twist is non-trivial, because 
of the additive mass renormalisation of the Wilson formulation, but achievable 
numerically in practise. 



3 The FNAL heavy quark formalism for Wilson 
fermions 

The FNAL formalism for heavy quarks was originally described by [T5]. Further 
developments of the FNAL lattice heavy quark formulation are described in 
1351 fT9l l20l [37] Here I review the free field calculation from El-Khadra et al. [T5] 
as warm up to adding a twisted mass term. 

The formalism starts with the quark propagator in momentum and time. 

at* -a -Et sinh ^7osign(t) - i%pi + m q + 1 - cosh E + \pipi 

S{t ^ = e 2Z^kE (9) 

where pi = sin(apj) and pi = 2sin(^ ! | i ). and Z2 is calculated to be 

Z 2 = 1 + m q a + ^p 2 a 2 (10) 

Equation is only valid for t > 0, because there is an additional term at t — 0. 
As am q gets very large, I/Z2 will get very small and this will cause problems with 
the dynamics. 

In the FNAL formulation, when the quark mass gets heavy, the dispersion 
relation of the heavy quark gets modified to 

E 2 = M? + ^±? + ... (11) 
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Mi = m q - -m 2 q + -m\ + O(mJ) (13) 



where M\ is known as the rest mass and M 2 is called the "kinetic" mass. Another 
closely related way to measure the deviations of the lattice dispersion relation 
from the continuum one is via the "speed of light" [38] . 

Mi = -ln(l +m q a) (12) 
Expanding the rest mass in terms of m q a gives 

1 2 1 
—m„ H — 5 

2 q 3 

The second term in the expansion is the leading b m improvement term in the AL- 
PHA formulation of the clover action [39]. The connection between the ALPHA 
formulation of the clover fermion action with the FNAL formulation of the heavy 
fermion action is demonstrated at one loop by Mertens et al. [10] for the quark 
mass. 

The M 2 kinetic mass is extracted using [TS] 

1 d 2 E 
M~ 2 ~ d Pl d Pl ' Pl=0 

For Wilson fermions the standard result is 

1 2 



(14) 



+ T-Z- (15) 



M 2 m q a{2 + m q a) 1 + m q a 
with an expansion in the quark mass: 

1 7 

M 2 = am q - -(am q ) 2 + (am q ) 3 - -(am q ) 4 + 0((am q ) 5 ) (16) 

The KLM factor is defined ]22j via equation [T71 

Z 2 KLM E<° I V^MO) | 0) latt = / rf 3 x(0 I V(^(0) | 0) cont (17) 

x J 

hence 



jklm = yZ 2 {p= 0) = y/1 + m q (18) 

There is a prescription for constructing amplitude factors for operators that are 
extended in space and time, such as conserved currents [HI [27] that I do not 
discuss. 



4 The FNAL heavy quark formalism for twisted 
mass QCD 

For the twisted mass formulation I first match the twisted heavy action onto a 
twisted continuum fermion action in equation HJ then rotate back to the standard 
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continuum action in the physical basis. This two step procedure seems more 
natural than trying to match the action in the twisted basis back to the continuum 
Dirac action in one step. I add the superscript T to show that the quantities are 
for the twisted action. 

The quark propagator in time and spatial momentum for Wilson twisted mass 
fermions has been written down by Cichy et al. [12] in the twisted basis, as part of 
their study of the pion and nucleon correlators in free field theory for the twisted 
mass action and a variety of actions that obeyed the Ginsparg- Wilson relation. 

S ^ ^ = 2Z T sinh ffr ^/fegnft) sinh ^74 - ij&i) 

+ [(1 - coshE) + am q + - iafx q j 5 r 3 )e~ E 1 (19) 

The rest mass (Mi) is obtained [321 [12] from the energy (E T ) at zero three 
momentum: 

„ a 2 m 2 n + a 2 u 2 

cosh M? = 1 + — ^ 20 

1 2(1 + am q ) v ; 

At maximal twist, equation |2"U1 shows the pole mass Ml is a function of the (an q ) 2 
so there will be no dependence on odd powers of the lattice spacing, consistent 
with the general symmetry arguments [12J. The corrections to the continuum Mi 
= afi q are much smaller than for the Wilson action in equation [13j 

Ml = afi q - ^(a^) 3 + g|j(a^) B + 0((a/i,) 7 ) (21) 



The kinetic mass (Mj) is 



_ ^A + (a N ) 2 
2 _ afMq 2 - (a//,) 2 1 ] 



The expansion of Mj in terms of \i q 
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Ml = afi q + -(a/i g ) 3 + Y^( a ^) 5 + °(( a ^) 7 ) ( 23 ) 



The mass dependent amplitude is 



zJ< LM = \/zI(p = o) = v/TT^; (24) 

At maximal twist Zl LM is 1, because m q is tuned to zero. Is is surprising that 
Zrlm IS independent of fi q , because I would have naively expected an expression 
that depended on the twisted mass fi q , but with no 0(a) errors. In the calculation 
Zrlm is independent of /x 9 , because as Shindler [12] notes, the \i q term and the 
Wilson terms "point in different directions" so don't interfere. However it would 
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be "cooler" to have a deeper more theoretical argument. The results for the Mi 
and M 2 masses do show a dependence on fia, beyond the continuum result, but 
with no 0(a) terms as expected. 

The KLM factor was originally obtained as part of deriving the transfer matrix 
for Wilson fermions [US EE] ■ The derivation of the transfer matrix was extended 
to twisted mass QCD by Frezzotti et al. [321 [T2"] . The normalisation of the fields in 
the derivation of the transfer matrix depended on a matrix called " B" in equation 
13 in [13]. The equvalent "B" matrix for twisted mass QCD is the same as for 
Wilson fermions and independent of the \x q mass [321 [12], and so is consistent 
with the KLM factor being independent of fi q at tree level. 

Although this analysis is focused towards twisted mass fermions, it is inter- 
esting to try and understand the \i q independence of the KLM factor. One way 
of getting some insight it to look at lattice actions with more symmetry such as 
those with Ginsparg- Wilson symmetry [HI 1451 146] , or those with a remnant of 
chiral symmetry such as improved staggered actions [47]. 

Liu and Dong have studied 0((amg) 2 ) and 0((a 2 mgAQ C£) ) effects in renor- 
malisation constants and the dispersion relation in numerical data [lElH^]. They 
found that the variant of the overlap action they used had lattice errors under 
control if they kept arriQ < 0.5. This numerical work suggests that a KLM factor 
for overlap fermions does depend on the heavy quark mass, although there are 
no O(amg) corrections as expected. The work by the TWQCD collaboration 
uses the overlap action with much larger masses to study mesons containing the 
bottom quark [5"U] . 

Aarts and Foley [51] have studied an overlap operator [431 146] in free field 
theory. They [5]] find an overall mass dependent renormalisation factor for the 
quark propagator that suggests a mass dependent KLM factor. However, there 
are a wide variety of different solutions to the Ginsparg- Wilson relation, some of 
which will have a different mass dependence. Liu and Dong [IB] discuss one choice 
that may have good properties in the heavy mass limit. It would be interesting 
to see if a KLM factor could parameterise the numerical data of [481 149] using 
methods in [5T], [52] . 

The tree level mass corrections to the improved staggered action called HISQ 
were considered by the HPQCD collaboration [47]. The coefficient of the Naik 
term was tuned at tree level to obtain a speed of light of one, up to errors of 
order 0((am q ) 12 ). The wave function renormalisation at tree level reported by 
HPQCD [47] for the HISQ action with the Naik term corrected with a mass 
dependent factor, had an explicit but weak dependence on the quark mass. 

From considering the HISQ and overlap actions above, it is unusual to have 
a wave function factor that does not depend on the physical quark mass. 
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5 Automatic 0(a) improvement and the FNAL 
formalism 



One very clever proof for the automatic 0(a) improvement [TT] . used that the 
action in the physical basis (Si att F m equation [25] with £=1 and r s =l) was in- 
variant under the symmetry V x T>d x (fi q — > — fi q ) where (xp = (—x,t)). The 
automatic 0(a) improvement of the twisted mass Wilson action, based on the 
above symmetry, has been tested numerically in quenched QCD [53j EH EH] and 
n/ = 2 unquenched QCD [T3] , 

The "P symmetry transformation is defined by 

U (x) -> C/o(arp) 
£4(x) — ► Ul(x P -ak) 

^P(x) -> 7o^(^p) 
V>(x) -> ^(x P )7o 

and the symmetry is defined by 

[^(x) -> C/J(-z - o/i) 
^(x) -> e 3i " /2 ?/>(-x) 
V>(x) -> e 3i " /2 ^(-x) 

The argument in [TT] only required that the action was invariant under the 
group PxD^x (jj, q — > — and not that the action is also invariant under the 
full hypercubic group. To connect with the proof in [TT], I consider the twisted 
mass action in equation [6] at maximal twist rotated into the physical basis in 
equation | 



Siat,F = V(aO(7oA) + (liDi + fi q - ir 3 7 5 (-^Do - r s -D^ + m cr )^(x) (25) 

X 

where m cr is m q tuned to the critical mass from setting the PCAC mass to zero. 

The action Si a t,F in equation EHl with arbitrary ( and r s parameters is also 
invariant under V x T>d x (fi q — > — fi q ), hence it should be automatically 0(a) 
improved. The clover term should have a coefficient that is an odd power of the 
quark mass. The two parameters: ( and r s , need to be tuned for the twisted 
version of the heavy quark action, but should only be an even power of the quark 
mass. The recent work on an improved Fermilab heavy quark action included 
dimension 7 operators |37J, so adding a twisted mass term, may help with the 
design of more highly improved heavy quark actions. 

I don't see any simple connection between the symmetry V x V d x (ji q — > —fi q ) 
and getting a mass independent KLM factor in equation | 
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6 Conclusions 



1 have discussed the inclusion of a twisted mass term with the FNAL heavy 
fermion action at tree level. This is useful for the analysis of existing Uf — 

2 twisted mass lattice QCD calculations with heavy masses, and heavy quark 
calculations using the Osterwalder-Seiler action [3U1 EI] on configurations with 
2+1+1 flavours of sea quarks. One surprising thing about the KLM factor for 
twisted mass QCD was that it was independent of the twisted mass at tree level. 
To estimate the order of magnitude of the various improvements terms I use the 
numerical values a s ~ 0.24 and arriQ ~ 0.26 for the f3 = 3.9 data set, with a lattice 
spacing of 0.0855 fm, from the ETM Collaboration |7|. For the twisted quark 
action with heavy mass ttlq, automatic 0(a) improvement means that the leading 
error should be 0((arriQ) 2 ) which is approximately 7%. The mass independence 
of the KLM factor for twisted mass fermions implies that the leading corrections 
are 0((a s amQ) 2 ) and numerically about 2%. The preliminary numerical results 
from lattice QCD calculations with heavy quark from twisted mas QCD seem to 
have larger errors than the above estimates |56| . 

I showed that the symmetry that protects the twisted Wilson action from 
0(a) corrections, should also protect an action where the hyper-cubic invariance 
is broken, as used in the FNAL heavy quark action. Some quantities such as 
the hyperfine spitting in charmonium are known to be sensitive to the value of 
the clover coefficient and lattice spacing errors [SZl EHl ESI [HO] so automatic 0(a) 
improvement should be useful. 

Since twisted mass QCD has no 0(a) errors, it is in principle possible to 
use lattice calculations at three different lattice spacings and take a consistent 
continuum limit for calculations that include heavy quarks. It may be useful 
to supplement the "brute force approach" with an estimate of systematic errors 
from the FNAL heavy quark formalism. 
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